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The problem of vortex structure in a single Josephson junction in an external magnetic field, in the 
' absence of transport currents, is reconsidered from a new mathematical point of view. In particular, 

we derive a complete set of exact analytical solutions representing all the stationary points (minima 
and saddle-points) of the relevant Gibbs free-energy functional. The type of these solutions is 
determined by explicit evaluation of the second variation of the Gibbs free-energy functional. The 
' stable (physical) solutions minimizing the Gibbs free-energy functional form an infinite set and 

are labelled by a topological number A'^„ = 0,1,2,.... Mathematically, they can be interpreted 
as nontrivial "vacuum" (A^„ = 0) and static topological solitons (A^„ = 1,2, . . .) of the sine-Gordon 
equation for the phase difference in a finite spatial interval: solutions of this kind were not considered 
in previous literature. Physically, they represent the Meissner state (A''„ = 0) and Josephson vortices 
(A''^ = 1,2,...). Major properties of the new physical solutions are thoroughly discussed. An exact, 
5-H ' closed-form analytical expression for the Gibbs free energy is derived and analyzed numerically. 

Unstable (saddle-point) solutions are also classified and discussed. 
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I. INTRODUCTION 

^ ' In this paper, we reconsider the old physical problemi^ of equilibrium vortex structure in a single Josephson 
junction from a new mathematical point of view. The necessity of such reconsideration is motivated by the fact that, 
in spite of significant contribution by a number of authors (see, e.g., Refs4i^iS*i^), the problem did not find an exact 
and complete analytical solution in previous theoretical literature. In particular, the basic question, namely, "What 
should be called an equilibrium Josephson vortex in precise mathematical terms?" , remained unanswered. Over years, 
this issue has been a source of considerable misunderstanding. For example, there still exists a wide-spread erroneous 
belief®'^- that Josephson vortices "do not form" in "small" junctions with W <C 2Aj, where W is the length of 
the insulating barrier and \,j is the Josephson lengthi^*^ (In reality, Josephson vortices do form for arbitrary small 
(-\^ W > 0, provided the externally applied magnetic field H is sufficiently high: see section V of this paper.) 

To clarify the situation, we consider the simplest case of a junction in a constant, homogeneous external magnetic 
, field, in the absence of externally applied currents. Relevant geometry in presented in Fig. 1. In particular, the x 
' axis is perpendicular to the insulating layer (the barrier); the y axis is along the barrier. The barrier length W = 2L 
is assumed to be arbitrary: < W < 00. A constant, homogeneous external magnetic field H is applied along the 
axis z. H = (0, 0,H > 0). Full homogeneity along the z axis is assumed. 

The difficulties of most previous approaches to the problem arose from the incompleteness of a traditional 
. formulationSiiii^ in terms of a static (time-independent) sine-Gordon equation 

o: ^''^ 1 • ^ m 



for the phase difference (p = cj){y) and physical boundary conditions 



c5 ■ — {±L) = 2edH {h = c=l). (2) 

dy 

Indeed, equations Q constitute an ill-posed boundary- value problem.^^ in the sense that they do not meet the 
requirement of the uniqueness of the solution. [This is a consequence of nonlinearity of Eq. Q and of the physical 
boundary conditions being imposed on ^ instead of cf).] Thus, although the general solution to was well-known^ 
the constants of integration specifying all physically observable configurations could not be unambiguously determined 
from boundary conditions Q alone. 

An obvious way to overcome the problem of incompleteness lies in an observation that Eqs. (Q), |(2J) are nothing but 
stationarity conditions of a relevant Gibbs free-energy functional! ^'^1^^ Therefore, among all solutions to ^ that satisfy 
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FIG. 1; The geometry of the problem: t is the thickness of the barrier; W = 2L is the length of the barrier; A is the London 
penetration depth; d = 2A + 1 is the width of the field-penetration region (shaded). The external magnetic field H is directed 
into the plane of the figure. 



© at a given H, one has to choose those that minimize (at least locally) the Gibbs free-energy functional. By general 
physical arguments, the stable (in a mathematical sense) solutions should represent an observable thermodynamically 
stable state (in the case of an absolute minimum), as well as observable thermodynamically metastable states (in 
the case of local minima) . Unfortunately, the issue of exact minimization of the Gibbs free-energy functional did not 
receive appropriate attention in the previous literature. (However, stability of some particular solutions was analyzed 
from a somewhat different point of view, e.g., in Refs^^^iSi.) 

The fact that a complete set of stable solutions to (Q) , © can be obtained in a closed analytical form has been first 
demonstrated in Refs^^i^^ within the framework of minimization of a certain wide class of energy functionals whose 
particular representative is the Gibbs free-energy functional of a single Josephson junction. Exact analytical solutions 
of Refsii^ii^ form an infinite set and are labelled by a topological (vortex) number Ny = 0, 1, 2, . . .. As is pointed out 
in RefsiiSii^, these solutions constitute a new class of particular solutions to the sine-Gordon equation that can be 
interpreted as nontrivial "vacuum" (for Ny = 0) and static topological solitons (for A'^„ = 1, 2, . . .) in a finite spatial 
interval. Solutions of this type were not considered in the previous literature on the sine-Gordon equation?i^iiSiii 

In what follows, we show that the new topological solutions resolve the problem formulated at the beginning of the 
Introduction. To make the paper self-contained and independent of Refsi^^ii^, we employ a new, more general method 
of the derivation of topological solutions. In particular, we start with a full set of solutions to ©, I©. In contrast to 
the approach of Refs.^'^'^'*, the selection of stable (physical) solutions is made by explicit evaluation of the sign of the 
second variation of the Gibbs free-energy functional, which allows us to establish general conditions of stability. 

In section II, we introduce the Gibbs free-energy functional and show that all its stationary points are either minima 
or saddle points. In section III, the second variation of the Gibbs free energy is discussed. We show that the sign of 
the second variation is related to the sign of the lowest eigenvalue of a certain Sturm-Liouville eigenvalue problem. 
In this way, we establish general conditions of stability. In section IV, a full set of solutions to ©, I© is obtained 
and classified with respect to stability. In section V, a full set of stable (physical) solutions to ©, ©is derived. 
These solutions comprise a Meissner solution {Ny — 0) and vortex solutions (iV„ = 1,2, . . .). Their major physical 
and mathematical properties are thoroughly discussed. An exact, closed-form expression for the Gibbs free energy is 
derived. This expression is analyzed numerically in the cases of a "large" junction and a "small" junction. In section 
VI, unstable (saddle-point) solutions to ©, © are classified and discussed. In section VII, we summarize the main 
results of the paper and make a few concluding remarks. Finally, in appendix A, we analyze stability in two singular 
cases (the Meissner solution in a semiinfinite interval and the single-soliton solution in an infinite interval). 



II. GIBBS FREE-ENERGY FUNCTIONAL: MAJOR PROPERTIES 



From now on, we will employ dimensionless units. Thus, the length scale along the y axis is normalized to the 
Josephson penetration length Aj. The magnetic field is normalized to the superheating field of the Meissner state in 
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a semi-infinite junction^ = (edXj) ^ . The energy scale is normalized to —j^- (In particular, the flux quantum 
in our dimensionless units is given hy = tt.) 

In terms of the dimensionless units, a phase-dependent part of the Gibbs free-energy functional per unit length 
along the z axis takes the formiii4 



' dy 



-L 



1 — cos (p (y) 





2' 


dy 





-2H [cl>{L)-^{-L)] 



Note that the last term in Eq. Q is, physically, proportional to the magnetic (Josephson) flux 



with 



$j = -[0(L)-,^(-L)], 



(3) 



(4) 



(5) 



being the corresponding local magnetic field (equilibrium or not) . 
The stationarity condition of the Gibbs free-energy functional l|2Jl, 



' dy' _ 



-0, 



reduces to the static sine-Gordon equation for the phase difference 



dy'^ 



(6) 



(7) 



and boundary conditions 



dy 



(±L) ^2H >0. 



(8) 



[Compare with Eqs. 1^ and respectively, where dimensional units are employed.] 

The main property of the functional Q follows from the fact that it belongs to the class of regular functionals, 
i.e., satisfies a necessary condition of the minimum;^ - hence, all stationary points of (|2J) are either minima of saddle 
points. [Unfortunately, the statement of Refii that all stationary points of |(3J) "are minima" is incorrect.] 

As already mentioned in the Introduction, boundary conditions (jSJ do not ensure the uniqueness of the solution to 
Eq. {Tj). Moreover, solutions to Q, ©, as stationary points of the Gibbs free-energy functional can correspond 
to either minima or saddle points of this functional. Taking into account that only stable phase configurations, 
corresponding to minima of Q, are physically observable, the problem one has to solve can be formulated as follows: 
Find all solutions to (O, © that minimize (at least locally) Q at given H > 0. To resolve this problem, we have to 
turn to sufficient conditions of the minimum^ which requires an analysis of the second variation of Q . 



III. SECOND VARIATION OF THE GIBBS FREE-ENERGY FUNCTIONAL AND A RELATED 

STURM-LIOUVILLE EIGENVALUE PROBLEM 

Let (j) = (j) {y) he a particular solution to (TJ, (|SJ) [i.e., a certain stationary point of (j3Jl]. The increment of the Gibbs 
free-energy functional Q in the vicinity of this solution, induced by variations 4> (y) — > (j) (y) -I- Scf) (y), where (50 has 
a continuous derivative that obeys boundary conditions 



7^(-i) = ^(i^) = o, 

dy dy 



(9) 



can be expanded in an infinite series 
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dy ' 



H 



dy 



k>3 



(10) 
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Here, 







L 






2' 




dy _ 


^ 1 dy 


cos <j){y) [5(1) {y)f + 


dSd) , . 





(11) 



and 



J dy COS 




-L 





'(y)] \ k>3. 



The sign of the second variation Hll|l in the expansion IjlOf) determines a type of the tested solution 
Three different cases are possible: (i) the case 



dS(j) 
dy 



> 



corresponds to a minimum of Q, (i.e., the solution is inside a stability region); (ii) the case 

>0, 



dy 



(12) 

Hy)- 

(13) 
(14) 



corresponds to a boundary of the stability region (a bifurcation point), IS when the solution can loose stability with 



respect to a certain variation; (iii) finally, if S^Q,g 



dy 



has no definite sign, the solution corresponds to a 



saddle point of ||2Jl, which means absolute instability. Therefore, in what follows, we concentrate ourselves on the 
evaluation of S^ftc- 

As is shown in variational theory of eigenvalues,^" the functional d^flc satisfies the following general relationiSi 



dy 



> Ho dy[S(j){y)f , 



where /io is the lowest eigenvalue of the Sturm-Liouville eigenvalue problem 



dy 



2 + cos 0(y )■(/' = M'/', y^[-L,L] 



f^(-X)^f^(L) = 0. 
dy dy 



(15) 

(16) 
(17) 



Equality on the right-hand side of relation (|15|l is achieved, when S(j) coincides (up to a factor) with the eigenfunction 
corresponding to /ig: i.e., dcj) = const ^pQ. As is clear from relation (|15|l . the three different types (i)-(iii) of the behavior 
of S'^V,G correspond, respectively to: (i) > 0; (u) — 0; (in) < 0. 

The properties of the self-adjoint operator, specified by the left-hand side of Eq. H16|l and boundary conditions l|17|) . 
are well-knownf2i In particular, its spectrum /i„ {n = 0, 1, 2, . . .) is discrete, real, infinite and bounded from below: 



/io < Ail < Ai2 < 



(18) 



(Note that in our case \fio\ < 1.) The corresponding eigenfunctions tjjn [n = 0, 1, 2, . . .) are real, mutually orthogonal 
and can be normalized: 



Li 



dyipn (y) V'm (y) = 0, n^m; < / dy [i/jn (y)] < ao 



(19) 



The eigenvalue number n determines the number of nodes of the eigenfunction ipn in the interval [—L, L] . Thus, the 
eigenfunction ipo can be considered strictly positive: 



V'o = -00 (y) > 0, ye[-L,L]. 



(20) 
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Significantly, the set of eigenfunctions {V'nlJ^o complete in the sense that any function / (y) that possesses con- 
tinuous derivatives up to the second order and obeys boundary conditions / (±L) = can be expanded in terms of 

Upon substitution of a solution of the boundary- value problem Q, © into (|16|l . this equation can be transformed 
into the well-known^'' Lame equation. Therefore, the boundary-value problem for ^ — can, in principle, be solved 
analytically. However, since we need to know only the sign of fxo, we will derive below a set of exact relations that 
will allow us to determine the sign of fiQ without explicit evaluation of /iq- 

First, we notice that the function 

X^f, (21) 
ay 

where (p = (j){y) is a particular solution to iQ, © (in what follows, we drop the bar over </)), satisfies the equation 



dy 



^J-f cos(/.(y)x = 0. (22) 



Combining l|21() . I|22() and (|16|l (with ij, = /iq, ~ tpo), using boundary conditions H17|l . we obtain the sought relations 

V'o (L) (£) - ^0 (-£) g i-L) (L)sin0(£)-^o(-L)sin(/)(-L) ^^^^ 
Mo = 1 = 1 • (2^j 

_/ d#o(2/)|f (y) / d#o(2/)i^(2/) 

[For definiteness, we choose the sign of V'o according to H2()(l .] In addition, we note that, if i/jq is known explicitly, /io 
can be found from the relation 

L 

J dyipo (y) cos (j) (y) 

f^o - . (24) 

/ dyipo (y) 

-L 

IV. COMPLETE SET OF SOLUTIONS AND THEIR CLASSIFICATION WITH RESPECT TO 

STABILITY 

We begin with the first integral of Eq. Q: 

1 ^ 

cos</> = C, -l<C<oo, (25) 



2 



dy 



where C is the constant of integration. Using 125|l . it is straightforward to derive the general solution to {TJ.^^ We 
write down this solution in the following explicit form: 
(I) -1 < C < 1: 

0± (y) = (2?^ -I- 1) ± 2 arcsin [k sn {y — yo, k)] , n = 0, ±1, . . . , (26) 



(II) 1 < C < oo: 



-K{k)<yQ<K{k), /c=i-t^, 0<fc<l; 



0± (y) = 7r(2n + l)± 2 am I ^—^,fc) , n = 0,±l,..., (27) 



-kK {k) < yo < kK {k) , k=-^—, Q<k<l. 

1 H- G 



6 



Here, the functions am u and sn u are the Jacobi eUiptic amplitude and the eUiptic sine, respectively; K (k) is the 
complete elliptic integral of the first kindiSi The choice of sign (i.e., 0+ or </>_) and allowed values of the constants of 
integration yo, k in (jSEJ, (|27|l should be determined from the requirement of compatibility with boundary conditions 
that we rewrite in a somewhat generalized form: 

^ (-L) = const > 0, ^ (L) = const > 0; (28) 
ay dy 

fi-L)-'^{L). (29) 
dy dy 

A. Solutions of type I 

Consider solutions of type I [Eqs. (|26|l ]. First, we note that, for this type of solutions, 

|0(L)-0(-L)| <2^ (30) 

for any Q < L < oo. Taking into account that 

d(j)± 



dy 



±2kcn {y-yo,k), (31) 



where cnu is the elliptic cosine, from H28|) we find that the appropriate solution is = 0+. If L < ^, the constant 
of integration k G [0,1]. In contrast, for L > ^, the allowed values for k are fc = and k e [fc„i,l], where km is 
determined from the condition 

K{km)^L. (32) 

Condition yields an equation for j/o, 

sn (L, k) dn {L, k) sn (yo, k) dn (yo, k) — 

(where dnu = -^amu), whose unique solution is yo — 0. Thus, solutions of type I, compatible with boundary 
conditions ||SJ), have the form 

(y) = TT (2n + 1) + 2 arcsin [k sn (y, fc)] , n — 0, ±1, . . . , (33) 

where k e [0, 1] for i < f , and fc = 0, fc G [km, 1] for i > |. 

An analysis of stability of solutions H33() is straightforward. Thus, for fc = we have 

<^(y) =^(2n+l), n = 0,±l,.... (34) 

The exact eigenfunction in H24|l is "ipo — const > 0, which immediately yields fio — —1. For fc ^ 0, the denominator in 
(1^ is positive, whereas 



^ (L) = -sn (L, fc) dn (L, fc) < 0, ^ (-L) = sn (L, fc) dn (L, fc) > 0. 
dy' dy' 



Therefore, fiQ < 0. In other words, solutions 133|) correspond to saddle points of © and, hence, are unstable and 
physically unobservable. 

B. Solutions of type II 

Consider solutions of type II [Eqs. (|27|l ]. In contrast to the solutions of type I [see H30|) ]. now 

[</)(£) -0(-i) I <oo (35) 

for any < L < oo. Using the derivative 

'^^^±ldn(y^,k]^ (36) 



dy fc \ fc 
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from H28|) we conclude that the appropriate solution is = with k G (0, 1]. Condition H29(l yields an equation for 
2/0, 

.„(f..)c„(f..).„(f,*)e,.(^,.).0, ,37, 

If 7^ pK (fc) {p = 1,2,.. .) in H37() . this equation has two different solutions: 

yo = -fci^(fc), yo = 0. (38) 
Correspondingly, solutions l|27(l . compatible with boundary conditions 0, split into two distinct sets: 

(l)eiy)^n{2n + l) + 2am(j-+K{k),k^ , n = 0,±l,..., (39) 



and 



0o(2/) = 7r(2n+ l) + 2 am (^|, A:) , n = 0,±l,.... (40) 

The meaning of the subscripts e (even) and o (odd) should be clear from the following: for solutions (|39() . we have 
(j)e (0) = 2n{n+ 1) (n = 0, ±1, . . .); in contrast, for solutions (jH, 0o (0) = tt (2?i + 1) (71 = 0, ±1, . . .). Using Eq. 
and its derivative. 



d3<^ , d(j) 

— = cos (j)— 
dy ay 



g - COS0— , (41) 



we find that the "local magnetic field" ^ at y — has a minimum for 0e and a maximum for 0o. 

It is interesting to note that the two sets of solutions H39|) . H4UI) are related by the Backlund transformations 



Id 1 ± . (j)oT' 

-— [0o ± (Pel = ; sm ■ 



2dy'" k 2 ' 

which, in the general case of a time-dependent sine-Gordon equation, is a hallmark of complete integrabilityii^*iSiii 
By virtue of the symmetry relations 

0e hy) = -0e (y) + 47r {n + 1) , 

0o(-y) = -0o (y) + 27r (2n + 1) , 

the eigenfunction -00 in (|lti|l . H17|) is necessarily symmetric with respect to reflection: ipQ {—y) — ipo (y). As a result, 
for 6(j) = const tpo, all the odd terms S^"^^^^Ig (m = 1, 2, . . .) in expansion H10() vanish by symmetry. 

Stability regions for solutions (|39|) and H4U|I can be established as follows. First, we note that the roots of the 
equations 

pkK{k)^L, p=l,2,..., (42) 
k = kp £ (0, 1], form an infinite decreasing sequence of bifurcation points: 

l>ki>k2> ... (43) 
Indeed, the second derivatives of (pe and 0o, respectively, read: 

d'^'l^e , . ^ n -ri;sn {l,k)cn (f,A:) d^(j)o , , ^ fV A (V 



Using H44|l . we find 



k . (44) 



(±L) = (45) 
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for k = kp {p = 1,2,.. .). Substituting (j) = (f>e,o with k = kp {p = 1,2, . . .) into (|23() . using (|45|l and the fact that 
the denominator in H23I) is strictly positive, we find fiQ = both for cj) = (j)^^ and (j) — (j)o- [As a matter of fact, the 
derivatives 



for k — kp {p ^ 1,2,...), coincide (up to a normalization factor) with the eigenfunctions ipOe,o 

subdivide the interval / = (0, 1] into an infinite set of semi-open 



dy ' ~ "-P 

corresponding to the eigenvalues fiQe.o — 0.] 
The bifurcation points k = kp (p = 1,2, 
intervals: 



I T 

^p=0^pj 



where 

/o = (fci,l]; Ip ^ {kp+i,kp], p^l,2. 
The index of these intervals p = 0, 1, 2, . . . can be expressed as 

4'e,o (L) - (t)e,o {-L) 



P 



2-E 



(46) 
(47) 

(48) 



where [. . .] stands for the integer part of the argument. As can be easily verified using (|23l) and (|44|l . /^oe > 
and ^oo < for fc e /q; ^oe > and ^oo < for fc € (™ = 1,2, . . .); /xpe < and /xqo > for A; e /2m+i 
(to — 0,1,2,...). [We again emphasize that the equalities /ioe,o — are realized only at the bifurcation points 
/c = /cp (p = 1, 2, . . .).] On these grounds, we conclude that stability regions for the solutions (f) = (pe are given by 
the intervals /2m (ti = 0, 1, 2, . . .), whereas stability regions for the solutions (f) — 4>o are given by the intervals /2m+i 
(to = 0,1,2,...). 

It is instructive to illustrate the above general results by explicit evaluation of /io = /io (k) for /c <C 1. In this limit, 
the lowest eigenvalue /io can be expanded in a power series of k: 



n>l 



where /Zq"'' (fc) is of order fc" (n = 1, 2, . . .). Here, we restrict ourselves to the evaluation of ^q''' (k): i.e., /iq ~ Mq^'' (fc) 
In view of the asymptotic expansions 



(1) 



,(1) 



w!, ^ ^ ^ o ^ , 1^ , 2y ky k'^ . 2y 

(j)e [y) « 27r(n + 1) + — - y - —sm—, n = 0,±1, .. ., 



(49) 



cj,,{y)^ni2n+l) + ^-^+'^8m'^, n = 0,±l,...., (50) 

it is sufficient to take a zeroth-order approximation to ■00 ^ V'o ~ V^o"'- -'^^ zeroth order in k, equations H16() . H17|) for 
■00 become 

^^=0, ^ (_L) = ^ (L) = 0, 

"2/ "y "2/ 

whose solution is 



^oL = const. (51) 



Upon substitution of l|I^ - (|FT|l into (|^ . we obtain 



Moe-^sm— , /i^j=-/csm— . (52) 

Expressions (|52|) immediately yield bifurcation points for k <^ 1: 

W 

kp = — , p=l,2,..., (53) 
np 

in full agreement with (|42|) . 

Summarizing, in this section we have derived a complete set of the solutions to |(7|) (both stable and unstable), 
compatible with boundary conditions (jSJ. The only stable solutions [i.e., corresponding to the minima of are 
those given by Eq. and Eq. I^UJl, where k G {'m = 0, 1, 2, . . .) and k € hm+i {m = 0, 17 2, . . .), respectively. 
In the next two sections both the stable and unstable solutions will be analyzed in more detail. 



9 



V. STABLE MEISSNER AND VORTEX (SOLITON) SOLUTIONS 



It is convenient to introduce a unified classification of the stable solutions, directly related to their physical inter- 
pretation. To this end, we introduce a new integer, the vortex (or topological) number = 0, 1, . . ., by means of the 
definition 



0(L)-0(-£) 
2tt 



I f , del) 

-L 



(54) 



[compare Eq. (|^ ]. and fix the so far arbitrary integer n = 0, ±1, ±2, ... in Eqs. and by the condition 

(/)(0) = TriVi,. (55) 

[Here, cj) — (j)e for Ny = 2m [m = 0, 1, 2, . . .), and (p = (j)o for — 2m +1 (to = 0, 1, 2, . . .).] Finally, using boundary 
conditions (jS)) for the determination of fc, we arrive at the desired form for the stable solutionsti^ii^ 



My) = TT (7V„ - 1) + 2 am (I + if (/c) , fc) , k = k (H) 



(56) 



dn[—,k\^ — kH — ' Ny = 2m (m = 0, 1, . . .) : 



(57) 



(/)o(2/) =7r7V„ + 2 am (|,fc) , k = k{H) 



(58) 



dn ( 7 ^ 



kH. 



7V„ = 2m+l (m = 0, 1,...), 



(59) 



where A'^ = corresponds to the vortex-free Meissner ("vacuum") solution, and = 1, 2 . . . correspond to vortex 
(soliton) solutions. The stability regions in terms of the field H take the form 



0<H<Ho, 7V„ = 0, 



(60) 



with 7Jjv„ implicitly determined by 



(iV„ + \)K 



H 



— Hn^L 



(61) 



(62) 



According to the results of section IV, we have S^ilc > within the whole semi-open interval (|60|l . Analogously, 
S^fto > inside the semi-open intervals (|61|l . whereas S^flo > at their boundaries. [Note that the upper bounds of 
the stability regions (|60|) . H61|) are also determined by the condition S^^Iq > 0.] 
Solutions H56|l - H62|) satisfy the symmetry relations 



i-y) - -0 {y) + 2^^. 



and obey the boundary conditions 



TT <(j> i-L) < 0, 2TrNy < (j> (L) < 27r [ Ny + 



1 



which ensures the fulfillment of the stability conditions 



d^(b , , d'^(j) , , 

ay dy^ 



(63) 



(64) 



(65) 
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As to major properties of the stable solutions, we want to emphasize the following. 

The stable solutions form an infinite set and their existence regions (|(jO|) . (|61|l overlap (at least for two neighboring 
states). This fact not only ensures that the stable solutions cover the whole field range < H < oo, but also proves 
that hysteresis is an intrinsic property of any Josephson junction, irrespectively of the value of < oo. However, 
overlapping is stronger for larger values of W and at 1 <C < (X) may involve several neighboring states. On the 
other hand, overlapping decreases with an increase of H. 

Soliton (vortex) solutions with iV„ = 1,2,... exist for arbitrary small < oo, provided the field iJ > is 
sufficiently high. We also note that the single-soliton (N^ = 1) solution appears at a finite (for any W < oo) field 
H = \/ Hq — 1 > , which should be contrasted with the case of the infinite interval (— oo. +oo): see Eq. ()68|l below. 

Solutions with iVu = 1,2,... are pure solitons only at iJ = Hf^ _^ ~ 1, when 4}{—L) — 0, 4>{L) — and 

5^^G > 0- In the rest of the stability regions (|61|l . when 5^VIg > 0, we have solitons "dressed" by the Meissner 
field. As a matter of fact, solitons (vortices) are confined to the spatial interval [— ^Z], where I is determined from 
the conditions cf) = 0, (/) = 2TTNy. In contrast, the intervals [—L, ~l) and {I, L] are "reserved" for the Meissner 
field. Nevertheless, although Josephson vortices always exist against a background of the Meissner field, solutions 
H56|l - Hfi2|l with Ny = 1, 2, ... by no means can be thought of as a mere superposition of the Meissner and the vortex 
fields, because the principle of superposition does not hold for the non-linear equation 

A solution with = 0, 1, 2, . . . solitons cannot be continuously transformed into the solution with A'^ + 1 solitons 
and vice versa. Therefore, any transitions between configurations with different vortex numbers A^^ are necessarily 
thermodynamic first-order phase transitions. 

Some examples of the stable solutions (for A^„ = 0, 1, 2), obtained by numerical evaluation of (|56|I - H62(I . are presented 
in Fig. 2. In several particular cases, known from the previous literature, solutions H5t)|) - (|62|l can be expressed in 
terms of elementary functions. For example, at iJ = and arbitrary L < oo, there exists only the trivial Meissner 
{Ny = 0) solution (po {y) = 0, y G [—L, L]. In the low-field limit < <C 1, the Meissner solution [equation with 
Ny — 0] reads: 

27? H 



By changing the variable y ^ y — L and proceeding to the limit L — > oo in equations H56|l . H57|l with A'^ = 0, we 
obtain the exact Meissner solution in a scmiinfinite intervali^ 

H exp ( — y) 

0ooo(y) = -4arctan-; , ?;e[0,+oo). (66) 

1 + V 1 — 

Equation immediately yields the upper bound of the existence of the Meissner state (or the superheating field of 
the Meissner state) in the semiinfinite junction: = Hg = I. The local magnetic field induced by (|65|l . 



^ ■ ' 2/e[0,+oo), (67) 

(1 + V 1 — 7i ^ j exp y — smh y 

vanishes at y = +00, as it should. 

By proceeding to the limit ±L ±00, — + 0, we find two different solutions: namely, the trivial Meissner (A'^ = 0) 
solution (y) = 0, y E (—00, +00), and the well-known single-soliton (A^^, ~ 1) solution in an infinite intervali^*iSiii 

<^ioo(y) = 4arctanexpy, y e (— 00, -l-oo). (68) 

The local magnetic field induced by (|68|l . 

h {y) — cosh^^y, y £ {—oo,+oo), (69) 

vanishes at y = ±00. Note that, from a point of view of the analysis of stability, solutions H66|) . (|68|l constitute 
singular cases: see appendix A. 

Of particular interest is the limit H ^ max{l,l/F}, which, physically, corresponds to negligibly small screening 
by Josephson currents. In this limit, the overlapping of states with different A^^, practically vanishes, and equations 
(EZjl-llSSJ become^ 

(-1)^" 

0e,o(y) - TTNy + 2Hy - ^-j^ [sin {2Hy) - 2Hy cos {HW)] , (70) 
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(a) 



(b) 



(c) 




FIG. 2: Stable solutions for L = 3.0. (a) The Meissner solution [N^ — 0): curves 1 and 2 correspond to the cases < H < Ho 
and H = Ho, respectively, (b) The single- vortex solution {Ny — 1): curves 1-3 correspond to the cases H = Hq — 1 , 
Hq — 1 < H < Hi and H = Hi, respectively; the vortex is confined to the spatial intervals denoted by vertical arrows, (c) 

The two- vortex solution (JV„ = 2): curves 1-3 correspond to the cases H = yjHl - 1, yjHl - 1 < H < H2 s.nA H = H2, 
respectively; the vortices are confined to the spatial intervals denoted by vertical arrows. 



where Ny = [-^^j , according to The distribution of the local magnetic field induced by (|7n|l is 

h{y) ^ H -^-^[cos{2Hy)- cos [HW)]. (71) 

Equations (TUJ, (|71|l exphcitly demonstrate the existence of solitons (or Josephson vortices) in the case W <^1. 

Upon substitution of H56|) . (|57|l into (|3Jl, we obtain exact, closed-form analytical expressions for the Gibbs free 
energy. It is convenient to write down these expressions in terms of the average energy density uj {H, k) = 



LU {H, k) = LOe {H, k) 5M„,2ra + '^o {H , k) 5N^.2m.+l (m = 0, 1, . . .) , 



(72) 
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k^k{H) 

{i-e)w 



8H 



W \ TT 

am I — + K (k) ,k 



2k 



Uo {H, k) = 2H'^ 



W 



8H (W , 



4fc2 
k^k{H), 



(73) 



where SN^,2m and (5iv„,2m+i are the Kronecker indices; E {u, k) and E (fc) are, respectively, the incomplete and complete 
elliptic integrals of the second kind.^** Note that in the limit H ^ max{l, W} [i.e., in the domain of validity of (|70() ]. 
equations (O, (|75|l reduce to 



I sin iHW)\ 
HW 



1 



8H^ 



coa^ (HW) 



(74) 



In Figs. 3(a) and 3(b), we present ui (H), obtained by numerical evaluation of (|72|l . (|73|l . for the cases of a "large" 
{L = 3) junction and a "small" (L = 0.3) junction, respectively. Figure 3(a) exhibits strong overlapping of neighboring 
stable configurations, whereas in Fig. 3(b) overlapping is practically invisible. The envelope of the energy curves for 
Ny = 0, 1, ... ,6 corresponds to the absolute minimum of the Gibbs free energy at a given H (a thermodynamically 
stable configuration). Parts of the energy curves that lie above the envelope in Fig. 3(a) correspond to local minima of 
the Gibbs free energy (thermodynamically metastable configurations). For better orientation in the physical situation, 
we have specified the upper bound of the existence of the Meissner state {H = Hq) and the first thermodynamic critical 
fieldiiSi^ {H = Hci). (The latter field is determined by the requirement that the Gibbs free energies of the states 
Ny = and Ny = 1 be equal to each other.) Both the fields, Hq and Hd, strongly depend on the length W: they 
increase (although at different rates) with a decrease of W. For example, for 1 <C < oo, they are approximately 
given by Hq Ri Hg = 1 and Hd ~ §, whereas for ^ 1 they practically coincide: iJ^i ~ ^^o ~ [see Fig. 3(b)]. 
Note that b y decrea sing the external field below H = Hd, one can still observe the single- vortex state down to the 
field H = ^ Hq — 1 < Hd [the abscissa of the left end of the energy curve for Ny = 1 in Fig. 3(a)]. 



VI. UNSTABLE (SADDLE-POINT) SOLUTIONS 

There are three different types of unstable (saddle-point) solutions to lO, ((SJ. Unstable solutions of the first and 
the second types obey symmetry relations 

^{~y)^-^{y) + 27rZ, (75) 
where Z is an integer. In contrast to the stable solutions of section V, this integer cannot be made to satisfy relation 

>(L)-0(-L)" 



Z 



2tt 



by any transformation (j) ^ (j) + 27rn [n = 0, ±1, . . .) and, thus, has no meaning of a topological number. In view 
of H75|l . the eigenfunction i/jq in H16|l . (|17|l is symmetric with respect to reflection: ip^ (— y) — tpa in)- However, in 
contrast to the stable solutions H5()|) - (|59|l . the first two types of unstable solutions are characterized by the property 

^W<0, '^i~L)>0, (76) 
dy^ ay' 

which results in /xq < [see (j^ ]. 

The first type of unstable solutions is represented by the set (|33|l . Without loss of generality, it is sufhcient to 
consider the case of n = 0: 

(f){y) ^ TT + 2 arcsin [ksn (y,k)], k^k{H): kcn{L,k) = H. (77) 

Solution H77|l exists in the field range < H < cosh^^L. At H = cosh^^i, it becomes 

0(?/) = 4arctanexp?/, ye[—L,L]. (78) 
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(a) " 




FIG. 3: The average Gibbs free-energy density ui (H) = for the cases of a "large" {L — 3.0) junction (a) and a "small" 

(L = 0.3) junction (b). Energy curves are shown for A''„ = 0, 1, . . . , 6. The upper bound of the existence of the Meissner state 
(-ffo) and the first thermodynamic critical field (Hd) are also shown. 



We emphasize that (|78|l is unstable for any L < oo, which should be contrasted with the stable single- vortex solution 
in an infinite interval, given by Eq. (|68|l . If L < ^, solution H77|) a,t H — degenerates into (p = tt. In contrast, if 
i > ^, aside from the solution ^ = tt, at = there exists a non-trivial solution 

= TT + 2arcsin[A:™ sn (y, fc™)] , (79) 

where km is implicitly determined by (|32|l . Solution H77|) is presented in Fig. 4(a). 

Unstable solutions of the second type can be obtained by prolonging solutions H5t)|) - (|59|) beyond their respective 
stability regions Ht)U|l . (|61|l . Several examples of such solutions are given in Fig. 4(b). 

Finally, unstable solutions of the third type do not obey (|75|l . Therefore, the eigenfunction i/jq in H16|l . (|17|l does 
not possess reflection symmetry: ipo (— y) ^ ipa (y). These solutions can be obtained from solutions (|39|l . (|40|l . taken 
at the upper boundaries of their respective intervals of stability, l2m with m = 1, 2, . . ., and l2m+i with to = 0, 1, 2, . . . 
(i.e., at the bifurcation points k = kp, p = 1,2, . . .), by a shift of the argument y ^ y — ak, where \a\ < K {k). This 
results in the properties 

^W-^(-^)>0, '^iL)^'^i-L)^0. (80) 
ay ay ay^ ay^ 

Relation ^ yields no < 0, because (L) < ipo {-L) for (±L) > 0, and Vo {L) > i^o {-L) for |^ (±i) < 0. 
Examples of unstable solutions of the third type are presented in Fig. 4(c). 



VII. SUMMARY AND CONCLUSIONS 



The main results of this paper can be summarized as follows. Mathematically, we have completely solved the 
ill-posed boundary value problem jT)), ijHJl and presented an exhaustive classification of the obtained solutions with 
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(a) 




FIG. 4: Unstable (saddle-point) solutions for L = 3.0. (a) The unstable solution of the first type: curves 1 and 2 correspond to 
the cases H = and H = cosh ~^L, respectively, (b) Unstable solutions of the second type: curves 1-3 correspond to different 
particular cases, (c) Unstable solutions of the third type: curves 1 and 2 correspond to different particular cases. 

respect to their stability. Our exact analytical treatment of the issue of stability can be easily extended to include 
more difficult cases, e.g., such as Eq. (UJ under more general boundary conditions^! or a coupled system of static 
sine-Gordon equationsiiii^ 

Physically, the complete set of stable solutions H56|l - (|62(l gives a clear final answer to the question what should be 
called the Meissner state and the vortex structure in Josephson junctions for arbitrary Vl^^ < oo and < iJ < oo. 
We have shown that the stability of physical solutions to Q , (jHJ is nothing but a consequence of soliton boundary 
conditions H64|l [or, equivalently, IjSSIl ]. In particular, the stability conditions IjtiSI) provide a simple, rigorous criterion to 
distinguish between physical (stable) and unphysical (unstable) solutions, which can be readily verified by comparing 
Fig. 2 with Fig. 4. 

To illustrate the difference between the results of our paper and those of previous publications, we turn to a 
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two-parameter expression 

0(j/) = ^ + 2am (^^^,fc^ . (81) 

Proposed in the previous literature2i2i^iSi2& as an "exact vortex solution", expressions of the type (|81|l by no means 
can be regarded as such in the absence of any expHcit definition of the constants of integration fc, yo and of stabiUty 
regions. As we have shown, aside from the Meissner solution [equations (|57|l . (|^ . with Ny — 0] and the 

actual vortex solutions [equations H5()l) - (|62|l with Ny = 1,2,...], the two-parameter family H81|l contains absolutely 
unstable, unobservable solutions of the second and the third types [see section VI and Figs. 4(b), 4(c)]. 

The reader should be warned against confusion between the single-soliton solution in the infinite interval 0i ^ [equa- 
tion (|68|l ] and its restriction onto a finite interval, the saddle-point solution H78|) [curve 2 in Fig. 4(a)]. Unfortunately, 
the unstable solution (|78|l is often erroneously interpreted in literature^" as a "Josephson vortex". 

Historically, the single-soliton solution oo served as the first example of vortex solutions in Josephson structuresi 
However, this solution cannot be realized experimentally in any realistic (with W < oo) physical system. In this 
regard, we want to emphasize once again that the actual physical single- vortex solution is given by Eqs. (|58|l . I|59|l . 
(EU and ^ with Ny = 1 [see also Fig. 2(b)]. 

Concerning saddle-point solutions, classified and discussed in section VI, they may play a certain role as channels 
of fluctuation-induced transitions from thermodynamically metastable states [local minima of (O] to a thermody- 
namically stable state [an absolute minimum of (Compare the decay of current-carrying states in narrow su- 

perconducting channels. ■^^) However, this issue asks for further investigation within the framework of a more general 
non-equilibrium approach. 

Finally, as is pointed out in Refs.^'^'^^ and in the Introduction, the exact analytical solutions (|56|l - (|62|l constitute a 
new class of static topological solutions to the sine-Gordon equations. We believe that they may find application not 
only in superconductivity, but also in other fields of modern nonlinear physics that deal with sine-Gordon equations. 
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APPENDIX A: ANALYSIS OF STABILITY IN TWO SINGULAR CASES 



In this appendix, we present an analysis of stability of the Meissner solution H6()|l and the single-soliton solution 
(|68|l . This analysis leads to singular Sturm-Liouville eigenvalue problems, formulated on the intervals [0, +oo) and 
(— oo, +oo), respectively. 



1. Meissner solution in a semiinfinite interval 



The Meissner solution (^ooo [equation H6()|l ] that exists in the field range < < 1 is a stationary point of the 
free-energy functional 



n 



Goo 



dy 



= / dy 



l-cos<j>{y) + ^ 



dcb{y) 



1 2 



dy 



2H(I){0). 



(Al) 



The second variation of (|A1|) at = 0o oo has the form 



d5<j) 
dy 



dy 



cos cj)ooo{y) [5(1) [y)f 



dy 



(A2) 



where the variation Scf) has continuous first derivatives and obeys boundary conditions 



, (0) = 0, <50(+oo) = O, ^(+(x)) = 0. 
dy dy 



(A3) 
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At iJ = 0, we have 0ooo = 0, and Eq. (jA2p immediately yields oo > for all allowed variations S(j). To 

evaluate the sign of S^flo oo in the interval < H < 1, we have to evaluate the sign of the lowest eigenvalue /ig of the 
problem 



cPip 

+ cos(j)ooo{y)ip = fj,tp, ye[Q,+oo), 



dy 



(0) = 0, ?/'(+oo) = 0, 



where the normalizable eigenfunction ipQ has no nodes in the interval [0, +oo) and can be considered positive. 
For fiQ, we have the following general relation [compare with (|^ ] 



Mo = 



V'o (0) 



(0) 



/ rf#o {y) ^ [y) 



d d> 

that holds in the whole interval < H < 1. The evaluation of j'^^'^ (0) yields: 



'Ooo 



dy' 



(0) = - 



2H 



1 + \/l -i72 

d^<t>o 



1 + Vl - H' - H' 



(A4) 
(A5) 



(A6) 



(A7) 



According to (0) < for < < 1, and ^^-^ (0) = for = 1. Taking into account that the 

denominator in ljA6|l is positive [we remind that = 2ft,], we conclude that /xq > for < < 1, and /io = for 

H =1. Summarizing the results for H = and < H < 1, we state: S'Ug oo > for Q < H < 1, and S'flc oo > for 
H = 1. 



2. Single-soliton solution in the infinite interval 

The single-soliton solution (j>i oo [equation (jSHJ] that exists at = is a stationary point of the free-energy functional 

1 



dy 



= / dy 



1 — cos (j){y) + 



\dcl^{y)] 


2' 


dy 





The standard analysis of S'^Q at 



bi oo requires evaluation of the lowest eigenvalue /Zq of the problem 



dy 



^ + cos 4>i^{y)->jj = fxip, ye(-oo,+oo), 



'0(±oo) = 0. 



As usual, the normalizable eigenfunction ipQ has no nodes in the interval (— cxd, -l-cxj). 

In this particular case, both and ipo can be determined explicitly. Indeed, the first derivative 



xiy) = — {y) 



y € (— oo, -|-oo) 



dy cosh y 

has no nodes and satisfies boundary conditions x (±00) — 0- Moreover, it obeys the equation 



dy 



— -f cos 01 00 (y) X = 0, y£{-(X),+(X)). 



Upon a comparison with l|A9|) . (|A10|I . we conclude that 

V'o {y) = const cosh^^y, /xq = 0. 



(A8) 



(A9) 
(AlO) 

(All) 

(A12) 
(A13) 
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Thus, S'^Q. > for (/> = (pioo, and the single-sohton solution H68|l corresponds to a bifurcation state. In view of 
translation symmetry of the problem, a shift of the argument y ^ y — yo, \yo\ < oo does not affect the stability of 
the single-soliton solution 0i oo , which should be contrasted with the situation for soliton solutions in the bifurcation 
state in the case of a finite interval [—L, L]: see the last paragraph of section VI and Fig. 4(c). 
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